Circuit theory for noise in incoherent normal— superconducting dot structures 
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We consider the current fluctuations in a mesoscopic circuit consisting of nodes connected by 
arbitrary connectors, in a setup with multiple normal or superconducting terminals. In the limit of 
weak superconducting proximity effect, simplified equations for the second-order cross-correlators 
can be derived from the general counting held theory, and the result coincides with the semiclassical 
principle of minimal correlations. We discuss the derivation of this result in a multi-dot case. 

PACS numbers: 74.50.+r 74.40.+k 73.23.-b 72.70.+m 
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Fluctuations of charge current in mesoscopic struc- 
tures are in general sensitive to the interactions and the 
fcrmionic nature of electrons. In multi-terminal setups, 
the geometry of the circuit is important for the cross- 
correlations, and in superconducting heterostructures, 
also the Andreev reflection, the superconducting prox- 
imity effect and transmission properties of NS interfaces 
need to be accounted for. 

The general theory for the full counting statistics of 
current fluctuations in multi-terminal structures was out- 
lined in Ref. The calculation of the second-order cor- 
relators using this theory can be simplified, from compli- 
cated 4x4 matrix equations to a Kirchoff-type system 
for scalar parameters, using an approach discussed also, 
for example, in Rcfs. In the incoherent case^*^ the 
result coincides with the semiclassical principle of mini- 
mal correlations^^ In this paper we show the derivation 
of this result in a multi-dot system, and consider a few 
special cases. 

The theory considers a network of normal (Tn) and su- 
perconducting (7s) terminals (T = Tn U 7g) and nodes 
(AT), connected by connectors. Each connector is 
described by its transmission eigenvalues ^ and each 
node j is characterized by a Keldysh Green function Gj , 
which is a 4 x 4 matrix in the Kcldysh(") ® Nambu(") 
space. In the quasiclassical approximation, assuming sta- 
tionary state and isotropicity, these are only functions of 
energy, G(e). 

The statistics of the current in the circuit is connected 
to the generating function S({xk\ker) of charge transfer, 
which can be found by solving transport equations for the 
Green functions. In the stationary case at zero frequency, 
the noise correlations Ski between the fluctuations 51k = 
Ik — (Ik) of currents flowing into the terminals k, I G T 
relate to it through^ 



where the Green function has the form 



Ski = 



, 2 t dxkdxi 



{Xj}=0 



(1) 

Here, to is the duration of the measurement, and the 
equality applies provided this is much larger than the 
correlation time of the fluctuations. 

The boundary conditions for transport are assumed 
such that the terminals are in an internal equilibrium, 



G 



eq 



R K 
6 A 



R = MT3 + VIT2 , K = Rh — hA , 
A = -n&ta , h = f L + T 3 f T . 



(2) 



Here, u = \e\ j^Je 1 — A 2 , v = sgn(e)\/u 2 — 1 are the co- 
herence factors, and A is the superconducting pair am- 
plitude. The functions /r(e) = 1 — fo(s) — fo( e ) and 
/i(e) = /o( — e) — /o(e) are the symmetric and antisym- 
metric parts of / (e) = [ e {^-^)/{k B T) + ^-^ where T ig 

the temperature and V the potential of the terminal. We 
assume V = in all S terminals to avoid time-dependent 
effects. For calculation of the statistics of the current, the 
counting held theory additionally specifies the rotation 

Gk(Xk) = e^/ 2 G k , cq e~ l * fK ' 2 , f K = f x ® f 3 , (3) 

at each terminal k, which connects the "counting fields" 
Xk to the Green functions. 

In circuit theory^ transport is modeled by the conser- 
vation of the matrix current at each node i 
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T n [Gj,Gi 



(4) 



The sum runs over all nodes and terminals (C = T U 
M): we assume the convention that = for i = j 
and disconnected points. This matrix is related to the 
observable charge and energy currents by 



i r°° 1 r°° 

(5) 

Their dependency on {xi}, in turn, describes the gener- 
ating function of charge transfer ti 



dS({ x *}) = --rEE Jjfc ({»}) d( 



iXk) 



(6) 



Determining the Green functions at the nodes from 
Eqs. (I.'SI Hi and finally applying Eqs. (|5I6|) . one can in prin- 
ciple find the distribution of the fluctuations in the cur- 
rent. However, the problem becomes considerably sim- 
pler if one is interested only in the second moment of this 
distribution, i.e., the current noise as given in Eq. 

We proceed calculating the noise by assuming that the 
superconducting proximity effect is negligible, so that the 
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anomalous parts (oc t\,T2) of the functions vanish in each 
node^i Then, one can expand the Green function at 
node j to the first order in the counting fields {xk}, in 
the Nambu-diagonal formiSiiS 



r _ / T3 2hjT 3 



keT 



bi 



(7) 

where b k {e) = ibl.(\e\), c? k = c? T + f 3 c^ and h = f L + 
Tsfr- This satisfies the quasiclassical normalization = 
1 up to second order in {xfc}- For the matrix currents, 
the above corresponds to the expansion 




J coh. 



of Eq. Q, where Iq, h,k an d I c ,k have the structure / = 
T3l(T3e), due to symmetries in the Nambu space. Here, 
Icoh.({Xk}) contains the off-diagonal Nambu-clcments, 
present if j corresponds to a superconducting termi- 
nal. In what follows, we neglect this coherent part of 
the current, assuming there are additional decoherence- 
inducing sink terms in Eq. This is valid provided 

that the Thouless energy describing the inverse time of 
flight through the node or the connector is much less 
than the characteristic energy scales of the problem, or, 
if there is a strong pair-breaking effect in the node, e.g., 
due to magnetic impurities. 

One can then consider expansion (jSJ) in detail, assum- 
ing a node i is connected to a node or terminal j. This 
yields four independent equations of conservation: 



0. 



b.k 



0. 



(9) 

in which It corresponds to the spectral charge current, 
II to the energy current, and the last two to a "noise" 
current, with the symmetric part defined as I c ,T,k(£) = 
Ic,k(s) ~\-Ic,k{~£)- The corresponding antisymmetric cur- 
rent I l c J L k is not needed, as we concentrate on the noise 
in the charge current. The spectral currents have the 
form 

4 = 9iM - It) , % k =9iM-bi), (10a) 

tf = /° , , ^jeT sa nd\e\<\A\, 
[9ij{f L -fh) otherwise. 

Thus, no energy current flows to the superconductors for 
lei < |A|. The fourth current is 
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C c,T, k ~ 9v{4-^)-{b\-h' k ){s ij {e)+s ij {-e)) , (ff) 

but it can be eliminated, see below. 

The factors gij and Sij (e) appearing in the expansion 
can be identified as the conductances and spectral noise 
densities characteristic of the connectors, and their exact 



form depends on whether the connector lies between two 
normal points (NN) or between a normal and a super- 
conducting point (NS). The expressions for the NS case 
are lengthy, so for simplicity we use here only the limits 
e <C A and e ^> A for superconducting Green's functions, 
effectively neglecting the exact form of the superconduct- 
ing density of states (DOS). In this approximation, for 
an NS connector at lei |A| or an NN connector, 
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(e) = -g" N {2-(r L + f T ) 2 -(f L + f T f 



T7*N N I fl 
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fT-fl-fr) 2 } 



(12a) 
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T«(l-T«). 



(12b) 



The result for an NS connector at |e| <C |A| is 



4 S & = L 2 4 S [1 - (fl)" ~ (/t) 2 + ^ S (/t) 2 ] , (13a) 



q NS 



E 



2(T^") 2 



(13b) 



FH S = (g NS 



r' e V 16 (^) 2 (1 T in (13c) 



as found through an expansion of Eq. (0J . Naturally, the 
results above agree with expressions for the noise gener- 
ated between two terminals, with being the differen- 
tial Fano factor&i£ 

The above equations are supplied with the boundary 
conditions 



4 = 0, f k (e) = f (e,V k ,T k ) . 



(14) 



where k and I are indices of terminals. These can be 
found by comparing expansion (7J to Eq. © (for N ter- 
minals), and by examining the expression for / (for S 
terminals). 

Finally, Eqs. lll'iKil) yield the result 



& = E 
= E 



de -I 



kj 



:r.i 



bi)(bi-(4) Sij (s), (15) 



for the correlations between terminals k and I. In the 
last step, we eliminated all c T k from the set of equations, 
which transforms the result to a sum over all connectors 
(i, j) in the circuit. 

The equations above have a simple physical interpreta- 
tion. The first two of Eqs. © describe the conservation 
of charge (T) and energy (L) currents at each energy in- 
terval [e,e + de}. With boundary conditions I|14ll2ll3fl . 
they yield distribution functions /£, f T of electrons at 
the nodes. In addition, one needs to solve from Eqs. @ 
the variable b l k , which characterizes the coupling between 
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terminal k and node i. It turns out that this quantity is in 
fact the characteristic potential introduced for semiclas- 
sical multiterminal calculations 4i Knowing /, the stan- 
dard two-terminal relations^ I|12I13|> gi ye the spectral 
noise densities in each connector, and Eq. I|15|) describes 
how these couple to the terminals. The final result is 
similar to the scmiclassical result in diffusive metals^i 
and coincides with the result in dot systems, see below. 

The assumption of all nodes being in the normal state 
resulted in a simple way to handle superconductors in one 
special case: first, it takes into account that no energy 
current enters superconductors at |e| < |A|, and second, 
assumes that other effects due to superconductivity are 
localized in only one connector, where both the conduc- 
tivity and the generated noise are modified. Our last 
approximation of a piecewise constant superconducting 
DOS simplifies the resulting expressions. 

We implicitly assumed above that there is no inelastic 
scattering which would drive the system towards equilib- 
rium. However, following Ref. Q, a strong relaxation of 
the distribution function in a node may be modeled by 
assuming that fj has the form of a Fermi function. In the 
case of relaxation due to strong electron-electron scatter- 
ing, the corresponding potential Vj and temperature Tj 
can be determined by taking the two first moments, J de 
and / dee of Eqs. ©: 



3e 2 



(16a) 



£ 9ij [T?-T* + C-\V?-Vf)] = Q. (16b) 

j&C\T s 

These describe the conservation of charge and energy cur- 
rents. If some of the nodes are in non-equilibrium, one 
can define the effective voltages and temperatures so that 
Eqs. (|16|) still apply for the whole circuit. In addition, 
one can model relaxation due to strong elcctron-phonon 
coupling by forcing Tj coincide with the lattice tempera- 
ture, so that only Vj need to be determined. 

It is illustrative to note that the quantum-mechanical 
counting-field theory agrees with the well-known princi- 
ple of minimal correlations, which is often used in semi- 
classical calculationsi^ In a typical model, one has the 
Langevin equations 



J2 i ij = o , /'•• = 9ij (v j -v i ) + sr* 



(17) 



where SP J are the microscopic fluctuations of the current, 
generated in the connector Eliminating voltages Vj 

at the nodes and assuming they do not fluctuate at the 
terminals, one finds the result 



6I k = J2(bi-V k )5P- 

(i,3) 



(18) 



for the fluctuations SIk in the current flowing to termi- 
nal k. Assuming SI lJ are independent and evaluating 
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FIG. 1: (Color online). Total Fano factor Ftot for a 
string of N + 1 nodes (inset), connected by M connec- 
tors with the transmission eigenvalue distribution p(T) = 
2g(n^T(l — T)) _1 of chaotic cavities. Results are shown for 
three types of relaxation in the nodes. Inset: M connectors 
in series. 



■| ({51k, SIi}), one finds Eq. I|15f) . This coincides with the 
prediction from the counting field theory, for an arbitrary 
circuit, provided it is understood that Sy should be eval- 
uated using the (average) distribution functions at the 
nodes. These may in general be in non-equilibrium, and 
should be obtained from a kinetic equation. Moreover, 
in the incoherent limit, the semiclassical result is correct 
also in the presence of superconducting terminals!^ 

The above discussion also clearly shows that an at- 
tempt to evaluate the higher correlators of noise using 
the principle of minimal correlations fails, as this cor- 
responds to truncating expansion J7J after the first two 
terms. The higher-order semiclassical corrections needed 
to fix this arc discussed for example in Ref. 

Consider now an example setup that consists of M — 1 
nodes between two terminals "0" and "M" (see inset of 
Fig. nj, and attempt to calculate its differential Fano 

factor, F to t = ^f-\i=o a t zcr0 temperature T = Tm = 
0. For simplicity of resulting expressions, we assume that 
all connectors are identical, sharing the same distribution 
{T£>} of the transmission eigenvalues. 

First, if both terminals are normal, the application of 
Eqs. H9I10I12I14I15|) is analogous to the semiclassical cal- 
culation presented in Ref. I id and yields the result 



3F-1 
3M 2 



for no relaxation, 
Fiot^NN = < jj , for e-ph relaxation, 

I ^ , for e-e relaxation, M — > oo. 

(19) 

This shows that the limit M — ► oo corresponds to the dif- 
fusive limit, due to the isotropicity of electron momentum 
assumed at the nodes. 

If terminal M is superconducting, and relaxation is 
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FIG. 2: (Color online). Cross correlation S13 between the 
two normal terminals changes sign as a function of contact 
transparency T, while 524 and S13 stay negative. Voltages 
are chosen Vi = V3 = V , V2 = Va = 0, and all connectors 
are assumed to be identical. Inset: A four-probe system. 
Terminals 1 and 3 are normal, 2 and 4 are superconducting. 
Connectors are assumed to have sub-gap resistances R, r and 
differential Fano factors F, f . 

negligible, we need to apply Eqs. (|9I10I12I13I14I15|) . In 

this example b 3 k are then straightforward to find, f 3 L = /£ 

and f T = f^ + fr^M for j = 1, . . . , M-l. Summation 
in l(T5|) then leads to a simple result 

_ 2 (F NS - ^)R 3 NS + ^jT^Fnn - ^)R 3 nn 
tot 3 + [-Rjvsr + {M — 1)Rnn] 5 

(20) 

Here Rnn, FnNj Rns and -P/vs are the resistances and 
differential Fano factors of the NN and NS connectors, 
as given in Eqs. (|12|) and (|13|> . The result applies also 
for M = 1, and in fact, for M = 2 it is valid even if both 
connectors have differing {T^}. In the limit M — > 00, 
the Fano factor again tends towards that of a diffusive 
contact, showing the doubling of the shot noised 



Similar calculation shows that for strong inelastic e- 
ph scattering one has -Ftot.c-ph ~ Fnn/M for large M. 
For relaxation due to e-e scattering, in turn, Eq. (|16b(l 
first gives the temperature profile Tj = [Tq + £ _1 (V^ 2 — 
V/)] 1 / 2 , where Vj = b 3 V . From Eqs. (jl5!12ajl one then 

finds that J^tot.e-e = for M — > 00, showing again the 

doubling of the noise. Numerical results for the behavior 
at smaller M are shown in Fig. ^ 

It is mostly straightforward to solve the current corre- 
lations in multite rminal N -S systems, also discussed for 
example in Refs. Il6ll7ll8l For the four-terminal setup 
shown in the inset of Fig. |21 one obtains^ 

§13 = -c ie (\Vx\ + \V 3 \ + \Vi + V 3 \) ~ c 2 e \Vi ~ V 3 \ (21) 
ci = 16(r V - R)4 [rR(l+F) + 2FR 2 + 2r 2 (l-/)] (22) 

02 = w+w [r(2r +R+ 2r2Rl) + 2r2/ (23) 

+ (2 + rR- 1 )(2r 2 + ArR + R 2 )F] . 

Here, R = R N n, r = Rns, F = F NN and / = F NS , 
and the result is valid provided \Vi\ -C |A|, i = 1,3, and 
T = 0. If the connectors are assumed diffusive (F = 1/3, 
/ = 2/3), Eq. J2D agrees with Ref. E One also finds 
that for / > 1, the cross-correlation (|21ll can be positive 
if R is small enough^ contrary to the case in normal- 
state circuits. For NS contacts with transparen cy T, t his 
is satisfied for < T < 2{y/2 - 1), as in Refs. Il6ll9i A 
different example, where all four contacts are identical so 
that R is not small, is shown in Fig. [21 

In conclusion, we discuss a simple model for the 
transmission of noise in multi-dot incoherent normal- 
superconducting structures, applying the microscopic 
counting field theory. The formalism produces the princi- 
ple of minimal correlations, and has strong analogies with 
the scmiclassical theory of noise in diffusive structures. 
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